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Previously derived expressions for the characteristic function of work performed on a quantum 
system by a classical external force are generalized to arbitrary initial states of the considered 
system and to Hamiltonians with degenerate spectra. In the particular case of microcanonical 
initial states explicit expressions for the characteristic function and the corresponding probability 
density of work are formulated. Their classical limit as well as their relations to the respective 
canonical expressions are discussed. A fluctuation theorem is derived that expresses the ratio of 
probabilities of work for a process and its time reversal to the ratio of densities of states of the 
microcanonical equilibrium systems with corresponding initial and final Hamiltonians. From this 
Crooks-type fluctuation theorem a relation between entropies of different systems can be derived 
which does not involve the time reversed process. This entropy- from- work theorem provides an 
experimentally accessible way to measure entropies. 

PACS numbers: 05.30.-d, 05.70.Ln, 05.40.-a 



I. INTRODUCTION 

During the last decade various fluctuation and work 
theorems [1-4] have been formulated and discussed. 
They provide information about the fully nonlinear re- 
sponse of a system under the action of a time depen- 
dent force [5-8] in contrast to linear response theories 
in which the response is expressed in terms of correla- 
tion functions of the unperturbed system [9-13]. More- 
over, these theorems were used to establish the second 
law of thermodynamics and to sharpen its formulation 
[14-16]. Fluctuation and work theorems have primar- 
ily been derived, numerically tested and experimentally 
confirmed for classical systems which initially, or perma- 
nently arc in contact with a heat bath [17-19]. The sys- 
tem then initially stays in a state described by a canon- 
ical (Maxwell-Boltzmann) distribution at the tempera- 
ture of the heat bath. Classical isolated systems which 
initially are in a microcanonical state were investigated 
in Ref. [20]. Quantum mechanical generalizations were 
proposed recently [21-26], but only for canonical initial 
states. The equivalence of canonical and microcanonical 
initial states was demonstrated for classical systems in 
the thermodynamic limit [27]. Little emphasis has been 
given to the impact of general initial conditions on the 
work performed on a quantum system by a classical force 
[28]. 

In the present paper we investigate the work performed 
by an external force acting on an otherwise isolated quan- 
tum system. The characteristic function of the work is 
shown to always assume the form of a correlation function 
of the exponentiated system Hamiltonians at the initial 
and final time, regardless of how the system was initially 
prepared. This characteristic function comprises all as- 
pects of the statistics of the work. In the special case of a 
microcanonical initial state the inverse Fourier transform 



leads to an expression for the density of performed work 
which directly leads to a quantum generalization of the 
microcanonical version of Crooks theorem. This theorem 
relates the probability densities of work of the processes 
and its time reversed partner process to the difference of 
entropies of the equilibrium states that correspond to the 
initial and final Hamiltonians. Such relations have been 
confirmed for classical systems by computer experiments 
[20]. In real experiments an active time reversal is not 
feasible. We yet present an "entropy-from-work" theorem 
that does not contain the time reversed process and still 
allows one to obtain the entropy of the equilibrium sys- 
tem for the system with the final Hamiltonian, provided 
the entropy of the system with the initial Hamiltonian is 
known. 

Apart from its use in numerical investigations such 
as molecular dynamics [29, 30] or microcanonical Monte 
Carlo simulations [31], the microcanonical ensemble is 
known to provide the valid description of isolated sys- 
tems in equilibrium [32]. It presents the proper statistical 
mechanical framework for the description of isolated sys- 
tems of finite size such as clusters of atoms [33] , atomic 
nuclei [34] or Bosc Einstein condensates [35] and even al- 
lows for phase transitions in finite systems [36-38]. 

The paper is organized as follows. In Sect. II the 
characteristic function of work is expressed as a corre- 
lation function of exponentiated Hamiltonians. The den- 
sity matrix which enters this expression is given by the 
initial density matrix of the system projected onto the 
diagonal elements with respect to the eigenbasis of the 
initial Hamiltonian. In the special case of a canonical 
density matrix the known form of the canonical charac- 
teristic function is recovered. In Sect. Ill we consider 
a microcanonical initial state, and derive microcanoni- 
cal quantum fluctuation and work theorems. Moreover, 
the probability densities and characteristic functions for 
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canonical initial states are shown to be related to the re- 
spective microcanonical quantities by properly weighted 
Laplace transforms. Conclusions arc presented with a 
final Sect. IV. 



II. CHARACTERISTIC FUNCTIONS OF WORK 



The response of a quantum system on the perturba- 
tion by a classical, external force can be characterized by 
the change of energy contained in the total system. The 
energy of the total system is determined by its Hamil- 
tonian H(t). It includes the external, time- varying force 
and therefore depends on time. We will consider the dy- 
namics of the system only within a finite window of time 
[£o,£/] during which the force is acting in a prescribed 
way, resulting in a protocol of Hamiltonians, which we 
denote by {H(t)} tf j a - A measurement of the Hamilto- 
nian H(t) at a time t will result in an energy which is an 
eigenvalue e k {t) of the Hamiltonian, i.e. 



sity matrix 



H(t)ip k ,\(t) = e k (jk)(p k ,\(t), 



(1) 



where A is a quantum number which labels a possible 
degeneracy of the eigenvalue efe(t). The eigenfunctions 
<Pk,\{t) can be chosen as normalized and pairwise or- 
thogonal. They span the eigenspacc of H(t) belonging 
to ek(t). The projection operator on this eigenspace be- 
comes 



,(*) = X)Ivm(*)Xvm(*)I- 



(2) 



Measuring the Hamiltonian at the respective times of 
measurement i and tf, one obtains as results eigenvalues 
e n (t ) and e m (tf) of the Hamiltonian H(t). The work w 
performed on the system is given by the difference of the 
measured energies, i.e. by 



w = e m (t f ) - e n (t ). 



(3) 



The energy values e n and e m arising from the measure- 
ments arc random quantities. Consequently, the observed 
value of work, w, too is a random quantity. The proba- 
bility p n , with which the particular eigenvalue e n (to) is 
observed in the first measurement, depends on the den- 
sity matrix p(to) which describes the state of the total 
system at time to . According to the laws of quantum me- 
chanics this probability is given by the expectation value 
of the projection operator P n (to) onto the subspace of 
cigenstates of the Hamiltonian H(to) with energy e„(to), 
i.e. by 



P,, 



TvP n (t )p{t ) 



(4) 



Immediately after this measurement the system is found 
in the corresponding state with properly normalized den- 



Pn 



Pn(to)p(t )P n (to) 
Pn 



(5) 



and then evolves in time according to 



Pn(t) = U t ,t p n U£ to , (6) 

where the unitary time evolution operator U tj ^ obeys 
the Schrodingcr equation 

iWU tit0 /dt = H(t)U tM , U t0tt0 = 1. (7) 

Note, that in general the time evolved states ip n ,\(t) = 
Ut,t <fn,\(to) are not eigenfunctions of H{t). An excep- 
tion are quasi-static changes of the Hamiltonian for which 
the adiabatic theorem holds [39] The second measure- 
ment, at time tf, produces an eigenvalue e m (tf) of H(tf) 
which occurs with probability 



p{m\n) = Tr P m (t f )p n (tf) 

= TrP m (t f )U tf , toPn U+ tto . 



(8) 



The probability density of work, Pt f ,t ( w ) can then be 
expressed as 

Pt f ,t (w) =y]6(w- [e m {tf) - e n (to)]) p(m\n)p n - (9) 



For the characteristic function of the work, 

GW«)= / dwe mw Ptf , t0 (w), (10) 



being the Fourier transform of the probability density of 
work, cf. e.g. [13, 40], we find 



E 



e iue m (t f ) e -iue n (to) 



where 



x TrP m {t f )U tf ,t PnU^ ito p„ 
= J2^e luH{tf) P m (t f ) 

xU tj , to e-™ H ^p n U+j p n 
= Tr j uH{tf) U tf , to e- iuH ^p{t )U+ M 
= Tre iuHlI ^ ) e- iuH( - to) P (to) 

= ^ e iuH(t f ) e -iuH(t )j 



H H (t f ) = U+ tto H(t f )U tfito 



(11) 



(12) 



denotes the Hamiltonian in the Hcisenbcrg picture. Fur- 
ther, we used the completeness relation J2 m Pm(tf) = 1 
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and introduced the density matrix 

P(to) = Z^PnPn = ^ P nP{to)Pn (13) 

with respect to which the average in the correlation func- 
tion (e 1 "^^' e~ luH ^ to ^) to is performed. This density ma- 
trix describes the state of the system projected onto the 
cigenbasis of the Hamiltonian at time t = to. It coin- 
cides with the initial density matrix p(io) only if p(to) 
commutes with the Hamiltonian H(to), i.e. 

p{t )=p(t ) <^> [p(t ),H(t )}=0. (14) 

In particular, if the system is initially in a canonical 
state the known form of the canonical characteristic func- 
tion results [25] . From this very result the Tasaki-Crooks 
fluctuation theorem and the Jarzynski work theorem fol- 
low readily [25, 26]. Notably, the characteristic function 
(11) assumes the form of a two-point quantum correla- 
tion as it is the case in linear response theory, despite the 
fact that it embraces the full nonlinear response on the 
perturbation by the applied force. 



III. FLUCTUATION THEOREMS FOR 
MICROCANONICAL INITIAL STATES 

A microcanonical initial state is characterized by an 
ensemble that consists of all eigenstates of H(to) with 
eigenvalues in a narrow band around an energy E. All 
contributing states have equal weights. In the ideal case 
of an arbitrarily narrow energy range around E the state 
of the system is given by the following density matrix 

PE (t ) = Lo E 1 (t Q )5(H(t )-E) 

= LU E 1 (t )J2 5 { e k(to)-E)P k (t ), ( 15 ) 

k 

where 

u E (t )=Tr5(H(t )-E) 

= E^Mi )-£) (16) 

k 

denotes the density of states of the system described by 
the Hamiltonian H(to), and gu = Tr Pk{to) the degener- 
acy of the energy ek{ta). Eq. (15) presents a rather for- 
mal expression for the density matrix. It can be regular- 
ized in the standard way by replacing the delta function 
by any smooth function 8 e (H(to) — E) which approaches 
the delta function in the limit e — > 0. The value of e may 
be chosen such that the width of S e (H(to) — E) is of the 
order of the accepted uncertainty of energy. It is also con- 
venient to take an approximating delta function with infi- 
nite support such as a Gaussian function in order to have 
a well defined density matrix for all energies E even if E 
falls in a gap of the spectrum. Since the width e is always 
assumed to be much smaller than all other relevant ener- 



gies in all following expressions the formal manipulations 
with delta functions also apply for the respective expres- 
sions containing the regularized delta functions. Com- 
bining eqs.(ll) and (15) we find the following expression 
for the characteristic function of work 

G tf ,t (E,u) =uj- 1 (t ) 

x Tre iuffH(t ^ er mHHi - to) 6{E - H(t )) 
= UJ E 1 (t ) 

xTre iuH H (t f )-E s(H(to)-E). (17) 

The inverse Fourier transform is readily performed to 
yield the probability density of work, 

Ptftt0 (E,w) = J ^e~ iuw G tftta (E,u) 

^^E^o) 

x Tr S(H H (t f ) -E-w) 8(H(t ) - E). 

(18) 

From this expression one may formally obtain the cor- 
responding classical result by replacing the Hamilton 
operators by the corresponding Hamilton functions and 
the trace by an integral over all possible initial states. 
Thereby, the Hamiltonian at the final time, Hn(tf) must 
be replaced by the Hamilton function depending on the 
initial phase space variables Zq via the solutions z(zo, tf) 
of Hamilton's equations of motion at the final time. This 
results in the expression 

pf fM {E,w)=u E \t Q ) 

x J dz S(H '(z(z , t /), t f ) - E-w) 

x S(H(z ,t ) - E) 

(19) 

which agrees with the result obtained on the basis of clas- 
sical statistical mechanics in Ref. [20]. In the classical 
formulation it is of course correct to replace the energy 
E in the first delta function by the Hamilton function at 
the initial time as dictated by the second delta function. 
This transformation is not possible in the quantum ex- 
pression confirming the recent observation that work is 
not a quantum mechanical observable [25]. 



A. Microcanonical quantum Crooks theorem 

The dependence on the work w can be shifted in eq. 
(18) from the first to the second delta function by intro- 
ducing the final energy Ef — E + w. For the product of 
the density of states we(*o) an d the probability density 
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Pt f t (-E 1 , w) one obtaii 



where 



Lo E (t )p tj ,to(E,w) = Tr 6 (H H (t f )-E f ) 

xS{H{t Q -E f +w)) 

= Tr 5 (H(t f )-E f ) 

x6(H H (t. )-E f +w) 

= 0JE+w(tf) Pt .t f (E + w, -w), 

(20) 

where, in going to the second line, we used the cyclic 
invariance of the trace and the unitarity of the time evo- 
lution in order to transfer the time dependence from the 
first to the second delta function. Here the operator 
Hiiito) denotes the Hamiltonian in the Hciscnberg pic- 
ture for the time evolution running in reversed order with 
a backward protocol from the final to the initial time, i.e.: 



H H (t ) = Ul tf H(t f )U to , tf 



(21) 



In the third line of eq. (20) the result was again written in 
terms of the probability density of work as a function of 
the initial energy E. This relation can be formulated as 
the microcanonical quantum version of Crooks' theorem 
[4, 20, 26] saying that 



Pt f ,t {E,w) _ e [s„(E+w,t f )-S u (E,to)]/k£ 

Pto,t f {E + w, -w) 



(22) 



It relates the ratio between the probability densities of 
work for a process and the corresponding time reversed 
process to the difference of the entropies which belong 
to microcanonical equilibrium systems with Hamiltoni- 
ans H(to) and H(tj). The entropies are defined in terms 
of the respective density of states via the standard sta- 
tistical mechanical relation 



S u (E,t) = k B lnwE(i), 



(23) 



where ks denotes the Boltzmann constant. Note that 
the similarly looking canonical Crooks theorem [4] con- 
tains a nonequilibrium entropy production a->oooks which 
is given by Tuj Clooks = TScmoks - Q = w — AF where 
Scrooks denotes the nonequilibrium entropy, Q the heat, 
AF the thermodynamic free energy and w the nonequi- 
librium work. This entropy production is distinctly dif- 
ferent from the statistical mechanical expression (23) for 
the microcanonical entropy entering eq. (22). 



Q{E,t) 



dE'uj(E',t) 



(25) 



denotes the number of states below E. For systems with 
short ranged interactions the two definitions are known 
to coincide in the thermodynamic limit, i.e. in the limit 
of infinite systems. For small systems though Sn(E) has 
been proved to be more advantageous as it is an increas- 
ing function of energy by definition. [29, 42]. 

We note that for an ensemble of initial states with en- 
ergies uniformly distributed up to the energy E, a fluc- 
tuation theorem of the form of eq. (22) can be derived, 
in which the difference of the entropies 5 W is replaced by 
the corresponding difference of Sn, i.e. 



Pt < / ,to(^.«0 



,[Sn(E+w,t f )-Sn(E,t )]/ki 



-w) 



(26) 



This follows along the same line of arguments leading to 
the relation (22) for the uniform initial density matrix 
p < (t Q ) which is given by 



p <E (t ) = n E (t )e(E - H(t )) 



(27) 



and which was introduced by Ruelle [43] in the context 
of the microcanonical ensemble. The probability of work 
pf s to (E, w) for this ensemble is obtained from the cor- 
responding microcanonical probability by an integration 
over the energy, i.e. 

p< tt0 (E,w)=n- 1 (E,t ) / dE'u;(E',to)pt f ,t (E',w). 

J -co 

(28) 

The fact that the microcanonical quantum Crooks the- 
orem (22) depends on the time reversed process seem- 
ingly restricts its practical usefulness. As opposed to 
computer experiments it is impossible to perform an ac- 
tive reversal of time, i.e. to let time run backwards, in 
real experiments. 



B. Entropy- from- work theorem 



Actually, a yet different statistical mechanical entropy 
can be defined. It is proportional to the logarithm of the 
number of states below the energy E rather than of the 
density of states at the energy E [41] . This entropy then 
takes the form 

S a (E,t) = k B \nn E (t) (24) 



The experimentally inaccessible probability density of 
the time reversed process though can be eliminated by 
first expressing the initial energy in terms of the final 
energy and the performed work, and next, by integrating 
eq. (20) over all possible values of the work. In this way, 
the density of states at the later time can be expressed 
by an integral of the initial density of states, weighted by 
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Pt f ,t (Ef ~ w,w), i.e. 

dw u) Ef -w(to)Pt f ,t (Ef -w,w) 



dw LO Ef (tf)Pt ,tf(E f ,-w) 

= <»E f (tf) (29) 

With the definition (23) the following relation between 
the entropy of the initial system and the unknown en- 
tropy of the final system can be established: 



dwe s »( B t-~M ) pt f ,t (Ef-v),w) = e s ^ E f^f\ (30) 



Note, that for a fixed final energy the weighting func- 
tion pt,t (Ef — w, w) is not a probability density of the 
performed work. In the following we demonstrate that 
the left hand side of eq. (29) can be written in terms of 
a properly defined average of the exponentiated entropy 
conditioned on the final energy. 

In general, from a single initial energy not all relevant 
final energies are likely to be reached, or may even be 
impossible to reach. Therefore, in an experiment the ini- 
tial energies have to be scanned over a sufficiently large 
range of values. For each initial energy a sufficient num- 
ber of experiments has to be performed in order that 
a reliable statistics of work can be compiled for conve- 
niently binned final energies. Based on such a statistics 
the probability density of initial energies E conditioned 
on the final energy Ef, p tfito (E\Ef), can be inferred. If 
the initial energies E are uniformly sampled in a range of 
energies of size En the joint probability of initial and fi- 
nal energies is given by p(E, E f ) = p tfA> (E, E f -E)/E R . 
According to Bayes theorem the conditional probability 
density p tft t (E\Ef) becomes 



Ptf,t a (E\Ef) — j dEp(E,Ej) 



p t/ , t0 (E,E f -E) 
JdEp tftt0 (E,E f -Ey 



(31) 



Consequently, the integral over the density of states can 
be formulated as an average over initial energies condi- 
tioned on the final energy, yielding 

(wE(to))E f = J dEu E (t )p tf , to {E\E f ) 

= N(E f r 1 

x J dwu Ef - w (t )p tf ,t (Ef -w,w) (32) 

where 

N(E f )= f dEp tftt0 (E,E f -E) (33) 
guarantees the normalization of the conditional aver- 



age. Using eq. (29) one finds the density of states 
aj Ef (tf) of a microcanonical system with the Hamilto- 
nian H(tf) represented by the average of the density 
of states uj E {to) with respect the conditional probabil- 
ity Pt f ,t (E\Ef) reading 



v Ef (t f ) = N(E f ){ue(to)) 



E f - 



(34) 



Expressing the density of states in terms of the entropy 
(23) one obtains 

e?~&ittV*B = N ( Ef )( e s„W(k B ) Er (35) 

We call eq. (35) the " entropy- from- work" theorem. In 
analogy to the Jarzynski relation it allows one to ex- 
tract equilibrium properties from non-equilibrium exper- 
iments, which are, in the case of the Jarzynski relation, 
the free energy and, in the present context, the entropy 
SUE)- 

Based on the relation (26) an entropy-from-work theo- 
rem for Sq (E) follows in close analogy to the correspond- 
ing relations (34) and (35) reading 



Q Ef (t f ) = N<(E f )(n E (t )) Ef . 



(36) 



and 



M E *>*tV'<° =N<(E f )(e s °W)< f/hB . (37) 



where the average ( - )^ is performed with respect to the 
conditional probability density 

pf fM {E\E } ) = N < (Ef)~ 1 pf }to (E, Ef - E). (38) 

and where 



N K (E f ) = / dEp< >to (E,E f -E). 



(39) 



C. Interrelations with the characteristic function of 
work for canonical initial states 



To establish a relation between the microcanonical and 
canonical work distributions we integrate the expression 
(17) for the micocanonical characteristic function of work 
weighted by the density of states and by a Boltzmann 
factor exp(— f3E) with inverse temperature (3. Performing 
this integration under the trace we obtain 



dEiOE^e-^Gt^t^E.u) 

= Tr e ™- f/ff( ^' ) e - iuH M e ~0H{t o ) 
= Z(t«,)Gf />t0 (u) (40) 

where 

Gf to (u) = Z fJ {t Q y 1 Tre mHH( - t f ) e- mHM e-l 3H ^ (41) 
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denotes the characteristic function of work for a process 
starting from a canonical density matrix [25] and 

Zp{t Q ) =Tre- fmM (42) 

the respective partition function. Hence, the canonical 
characteristic function of work is related to the micro- 
canonical one by a Laplace transform. The reverse rela- 
tion is readily obtained by the following inverse Lapace 
transform 

G t/ , t0 (E,u) -^(to)- 1 j c ^-Z p {t Q )ee E G^ {u) (43) 

where C is an inverse Laplace contour in the complex f3 
plane from — ioo + c to ioo + c. The real constant c must 
be chosen such that all singularities of the integrand lie 
to its left side. The density of states is related to the 
partion function by the standard relation 

Mk) = I Y~. Zpih) e PE - (44) 

Accordingly, the probability densities of work for mi- 
crocanonical and canonical initial states are also related 
by a Laplace transform, yielding 

p? f!to (w)=Z /3 (to)- 1 J dEuj E {t Q )e-P E PtfM {E,w) 

(45) 

p tfM (E,w)=u E (t )- 1 [ ^-Z (t o )e^^ to ( W ) 
J c 

(46) 

Finally we note that the Jarzinski relation can be ob- 
tained from the microcanonical Crooks theorem. For this 
purpose one considers the left hand side of the first line 
and the right hand side of the third line of eq. (20), mul- 
tiplies both sides with exp (— (3(E + w)) and integrates 
over all values of E and w. For the left hand side one 
then finds: 

J dE J dwer^ E+w ^ E (t )p tf ^(E,w) 

= J dwe-^Zp(t )pP fito (w) 
= Z p (to){e-^) (47) 
while the right hand side yields 

J dE J dw e-^ E+w ^ E+w (t f ) p to ,t f (E + w,-w) 

x J dwpt ,t t {Ef,-w) 
= Zp{t f ) (48) 



where we substituted the integration variable E by Ef — 
E + w. A comparison of the last two equations immedi- 
atly yields the Jarzynski relation [3, 25]. 



IV. CONCLUSIONS 

The expression for the characteristic function of work 
performed on an isolated quantum system by an exter- 
nal force was generalized for arbitrary initial states. The 
general structure of the characteristic function is given 
by a correlation function of the exponentiated system 
Hamiltonians at the first and the second measurement 
time. The quantum expectation and ensemble average 
are jointly taken with respect to a density matrix which 
results from the actual initial state immediately before 
the first measurement by means of a state reduction 
with respect to the energy eigenbasis of the then mea- 
sured Hamiltonian. Initial states which commute with 
the Hamiltonian at the initial time consequently are not 
modified. This form of the characteristic function holds 
irrespectively of a possible degeneracy of the spectrum of 
the system's Hamiltonian. 

For a microcanonical initial state expressions for the 
characteristic function of work and the respective prob- 
ability density were established. In the classical limit 
known expressions were recovered and the validity of the 
microcanonical Crooks fluctuation theorem was demon- 
strated for quantum systems. Moreover we formulated 
an entropy-from-work theorem which allows one to infer 
the unknown entropy of a system from a reference sys- 
tem with known entropy by means of a nonequilibrium 
experiment. In such an experiment the initial system 
with known entropy is deformed into the final system 
within finite time according to a prescibed protocol of 
Hamiltonians. In this context we want to emphasize that 
the entropy following from this theorem is based on ei- 
ther of the statistical mechanical definitions (23), or (24) 
which agree with each other for sufficiently large systems 
with short range interactions. In contrast, the canon- 
ical Crooks theorem contains a nonequilibrium entropy 
production that emerges from a relation involving free 
energy, heat and nonequilibrium work. 

We further note that the microcanonical distribution 
provides the appropriate description of an isolated system 
if no further information about its state is available even 
if the energy of the considered system is a priori unknown. 
By registering the result of the first energy measurement 
this ignorance is removed and the available information 
about the initial state is expressed without any bias by 
the microcanonical density matrix corresponding to the 
measured energy. Though, if the unforced dynamics does 
not only leave invariant the energy but also other quanti- 
ties such as linear or angular total momentum, then the 
adequate constrained microcanonical ensemble has to be 
considered as the proper initial state. 

Finally, by means of properly weighted Laplace trans- 
forms, relations between microcanonical and canonical 
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characteristic functions and probability densities of work 
were established. 
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